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We extend our subtractive-renormalization method in order to evaluate the 1 So and 3 Si- 3 Di NN 
scattering phase shifts up to next-to-next-to-leading order (NNLO) in chiral effective theory. We 
show that, if energy-dependent contact terms are employed in the NN potential, the 1 So phase 
shift can be obtained by carrying out two subtractions on the Lippmann-Schwinger equation. These 
subtractions use knowledge of the the scattering length and the 1 Sq phase shift at a specific energy to 
eliminate the low-energy constants in the contact interaction from the scattering equation. For the 

OJ=l coupled channel, a similar renormalization can be achieved by three subtractions that employ 
knowledge of the 3 Si scattering length, the 3 Si phase shift at a specific energy and the 3 Si- 3 Di 
\Q ■ generalized scattering length. In both channels a similar method can be applied to a potential 

with momentum-dependent contact terms, except that in that case one of the subtractions must be 
replaced by a fit to one piece of experimental data. 

This method allows the use of arbitrarily high cutoffs in the Lippmann-Schwinger equation. We 
examine the NNLO S-wave phase shifts for cutoffs as large as 19 GeV and show that the presence of 
linear energy dependence in the NN potential creates spurious poles in the scattering amplitude. In 
consequence the results are in conflict with empirical data over appreciable portions of the considered 
■ cutoff range. We also identify problems with the use of cutoffs greater than 1 GeV when momentum- 

dependent contact interactions are employed. These problems are ameliorated, but not eliminated, 
by the use of spectral-function regularization for the two-pion exchange part of the NN potential 
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I. INTRODUCTION 



Chiral perturbation theory (%PT) has now been applied to the problem of nuclcon-nucleon (NN) interactions for 
almost two decades. Since Refs. [l|, 0], chiral potentials have been deduced from xPT Lagrangians [H, 0, [H, S 0, II] 
through the computation of all NN-irrcduciblc diagrams that occur up to a given order in xPT. The resulting potential 
includes a variety of contact interactions, which parameterize our ignorance of physics coming from energies higher 
than the chiral-symmetry-breaking scale, A x . Some of these contact interactions occur in the 7rN Lagrangian, and the 
low-energy constants (LECs) that multiply these operators can be determined independently from the pion-nucleon 
scattering data [§, [l(| ■ More of a challenge — both conceptual and practical — is posed by the string of operators that 
represent the short-distance physics in the NN sector. These operators rcnormalize the divergent loops found when 
calculating the xPT potential, and so encode the contribution to observables from high-energy NN states. Each such 
contact interaction is associated with an unknown LEC. xPT power counting applied to the NN potential mandates 
that at a fixed order in the chiral expansion only a finite number of LECs enter the potential. For instance, at 0(P 3 ) 
in the NN S-waves (which will be where the bulk of our attention is in this work) there are five contact operators that 
must be considered: two act in the -"-So channel, two in the 3 Si channel, and one causes coupling between the 3 Si and 
3 Di partial waves. 
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For low partial waves the xPT potential is non-perturbative, and the standard approach involves use of the 
Lippmann-Schwinger equation (LSE) to reconstruct the full NN amplitude from the NN-irreducible part [ill, Ell 
US 111- ( c - f - Refs - E EE E3, El, EE HE E3, HI HE HI where parts of the potential are treated in perturbation 
theory.) Since the NN interactions obtained in xPT do not fall off as the nucleon momenta p go to infinity, a cutoff, 
denoted here by A, must be imposed on the intermediate states in the LSE. The LECs are then fitted to data for a 
variety of cutoffs, and should absorb any strong (power-law with power > or logarithmic) dependence of observables 
on A. A corollary is that the predictions of the effective theory should not depend on the particular quantity chosen to 
fix the LECs, as long as the kinematic point is within the domain of validity of the chiral effective theory. A potential 
which does not have these desirable properties (approximate cutoff and renormalization-point independence) is not 
"properly renormalized" . If the set of contact operators employed cannot achieve this we conclude that this "chiral 
effective theory" (xET) is unable to give reliable results, since its predictions are too sensitive to the treatment of the 
unknown short-distance physics at scale A. In Refs. [!, IE IE B 11] the potential V was computed to a fixed order, and 
then the NN LECs that appear in V were fitted to NN data for a range of cutoffs between 500 and 800 MeV. The 
resulting predictions — especially the ones obtained with the 0(P 4 ) potential derived in Refs. @,@] — are approximately 
A— independent and describe NN data with considerable accuracy. 

However, the range of A considered in these papers is quite narrow. Does the fact that these analyses only consider 
A < 800 MeV represent an intrinsic limitation on the xET approach? 

In the case of leading-order potential (which is 0(P°) and whose long-range part consists only of one-pion exchange) 
the answer has been shown to be "No" in the 1 So and and 3 Si— 3 D! channels. The problem can be properly renormal- 
ized for arbitrarily large A's with onl y co nstant contact interactions present in V in each channel — in accordance with 
the xPT power counting 0, HI, HEM HI, HI, [HE HE H3] x - However, once A > 600 MeV the xPT counting for the 
NN potential fails in waves with L > where the tensor force from one-pion exchange is attractive [13, H3, l3l[ [34l. [35| . 
That counting predicts no contact interaction in these channels, but a contact interaction is necessary to stabilize the 
LSE predictions at these cutoffs. Furthermore, in these partial waves it has been found that even when higher-order 
XPT potentials arc considered and the phase shifts obtained from the LSE are approximately cutoff independent, the 
results can still show significant renormalization-point dependence [35j . 

At next-to-leading order (NLO, 0{P 2 )) and next-to-next-to-leading order (NNLO, 0(P 3 )), the xPT potential 
includes diagrams with two-pion-exchange (TPE). Since the loop integrals there diverge as the square or cube of the 
momentum transfer |q| = |p' — p|, the TPE must be associated with contact terms up to order P 2 — just as one would 
have predicted using naive-dimensional analysis in powers of NN momentum. 

In Refs. [U, HH, H3, [HI the role of various different Q(P 2 ) contact terms in the LSE was examined and the rcnor- 
malization of the NN S-waves was discussed. Refs. [!, HU, H3] found that momentum-dependent contact terms had 
difficulty absorbing the cutoff dependence at larger cutoffs, and also that the contribution from the (supposedly 
higher-order) TPE was larger than one-pion exchange. (See Ref. [38| for a similar conclusion in the 1 So channel using 
a co-ordinate-space analysis.) Ref. [36| also considered an energy-dependent 0(P 2 ) contact interaction in the 1 So 
channel and showed that the use of a contact interaction linear in energy (which formally should be equivalent to the 
momentum-dependent interaction at the orders in the chiral expansion being considered here) leads to the appearance 
of resonances there for A — + oo. In cither case, the results of these studies suggest that for cutoffs larger than about 
1 GeV the contribution from contact terms dominates over that from the long-range part of the potential, raising 
the question of whether xPT is acting as a phenomenological fit form, rather than an effective theory. As argued in 
Refs. [T3, l35l. [38| . the chiral-symmetry-breaking scale could be responsible for this situation. 

This suggests that there may be a critical cutoff ~ A x above which it does not make sense to iterate the TPE inside 
the LSE 2 . It is the purpose of this paper to critically examine the behavior of the NN phase shifts predicted by x^T 
at NLO and NNLO in the % and 3 Si~ 3 Di channels over a range of cutoffs from A = 500 MeV to A > 2 GeV. In 
this way we hope to identify this critical cutoff, if it exists, and discuss the signatures of, and mechanisms responsible 
for, its appearance. 

In order to address these issues we extend our previously developed subtractive-renormalization method [32l . [35l | 
to evaluate the 1 So and 3 Si— 3 Di phase shifts with the higher-order xET potentials. If energy-dependent contact 
interactions are employed in these potentials we show how to relate the xET phase shifts to on-shell quantities. 
For momentum-dependent contact interactions we formulate the scattering equations in terms of two (one) on-shell 
quantities in the triplet (singlet) channel and then fit the single remaining unknown parameter to experimental 
data. This allows us to easily go to high cutoffs in the LSE and provides clean information about the effect of the 



1 In the 1 So channel the constant must be rriq-dependent [33l | in order to absorb a logarithmic divergence. This violates the power 
counting, but is of no practical consequence as regards the description of laboratory NN data. 

2 See [40| for a recent analytic approach which yields a similar conclusion. 
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renormalization point. Our method works for any long-range potential. Here we consider long-range parts evaluated 
using both dimensional regularization (DR) and spectral-function-regularization (SFR)[37|| at 0(P 2 ) and 0(P 3 ). The 
forms of these potentials were given in Ref. [35]. Since issues like differing extractions of the 7rN LECs ci, C3 and 
C4 could produce sizeable changes in the phase shifts we are not overly focused on the quality of our fit. Instead, 
our purpose is to see whether these potentials behave self-consistently with respect to cutoffs in the LSE. Thus, our 
strategy is to evaluate in detail the cutoff-dependence of the S-wave phase shifts, and examine the renormalization- 
point dependence. These features depend much less on the precise choice made for c±, C3 and C4. 

The structure of our work is as follows. First, in Sec. UH we introduce our subtractive-rcnormalization method 
and explain how to deal with both momentum- and energy-dependent contact terms in the singlet channel. We 
take as input to the method the 1 So scattering length, a s , and the phase shift at a particular energy, E* . Then, 
in Sec. IIII1 we extend our subtractive-renormalization method to the coupled-channels case and deal with constant, 
energy-dependent and momentum-dependent contact terms. While the constant contact term is, as in Ref. [32| . 
solved by one subtraction with at as input, for the energy-dependent contact term we pin down the three unknown 
constants with information with at, the generalized scattering length 020 that governs the low-energy behavior of ei, 
and the phase shift at a particular energy as inputs. In Sec. IIV| we discuss our results in the 1 Sq channel and their 
implications. In Sec. [Vj we discuss our results in the 3 Si- 3 Di channel and their implications. We summarize our 
findings in Sec. IVII 

II. SUBTRACTIVE RENORMALIZATION IN THE SINGLET CHANNEL 

The short-range part of the potential, i.e., the contact terms, is a parameterization of unknown high-energy physics. 
Before introducing our subtractive renormalization scheme, we need to decide what types of contact term we should 
consider. Throughout this work we will consider chiral two-pion-exchange 0, H, [H, [39[ as the long-range part of 
the NN potential. This is associated with contact terms up to 0(P 2 ). However, since for non-relativistic particles 
we have p 2 /M ~ E there is the possibility to consider energy-dependent, instead of momentum-dependent contact 
interactions. In particular, this may have certain advantages in terms of evading theorems that limit the impact of 
short-distance potentials on phase shifts [HI, HI]- Thus, in the 1 So channel we will consider the following three types 
of contact terms: 

(A) V SR,Q = A 

(B) v SR , (E) = \ + 1 E 

(C) V S R,0(P',P) = * + C2 (P 2 +P' 2 ), 

where A, 7, C2 are unknown constants. (Of course, the numerical value of A is different for each case.) The overall 
potential is then given by 

v (p',p;E) = v S R,o(p',p;E) +v LRt0 (p',p;E), (2.1) 

with the subscripts sfl an d lr standing for the short- and long-range parts of the potential. Specifically, in the 
channel we have 

v L r,o = (000\V C + W C - 3[V S + W S + q 2 {V T + Wr)]|000>, (2.2) 

with expressions for the S-wave projected central, spin, and tensor isoscalar and isovector potentials taken from 
Ref. [1, [!3] and collected in Ref. [35j| . 

A. The constant contact interaction 

The constant contact term can be evaluated by one subtraction with the 1 So scattering length a s as input [32| . 
This is possible because knowledge of the on-shell amplitude plus the long-range potential yields the fully off-shell 
amplitude, as long as only a constant contact interaction is present [32l. l43l. [44| . The step from the on-shcll amplitude 
at an arbitrary energy E* , tip^p^E*) (with p$ = ME*), to the half-off-shcll amplitude at the same energy is 
explained in Appendix [ATI The symmetry property of the amplitude t{p,p^;E*) = t(pQ,p; E*) then means that the 
same sequence of manipulations can be used to construct t(p' ,p; E*). With the fully-off-shell amplitude at a single 
energy in hand, it is well known how to construct the i-matrix at any arbitrary energy [26| . In operator form we have 



t(E) = t(E*) + t(E*)[g (E) - g Q (E*)]t(E), 



(2.3) 
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where go is the free resolvent of the Lippmann-Schwingcr equation, 



(Here M is the nucleon mass and E + = E + is, with s a positive infinitesimal. 



B. Energy-dependent contact interaction 

In this case there are two unknown constants. The strategy is to use one subtraction (as in Ref. [13]) to first 
eliminate the constant part of the contact term. For the energy-dependent contact term we then perform a second 
subtraction to eliminate the remaining unknown. 

We start from the partial-wave LS equation, which is written explicitly in this channel as 

p , , , ^2 f A dp"p" 2 Vo(p',p";E)t (p",p;E) 
t (p,p;E)=v (p,p;E) + -Mj o pg + fe-;/" ■ ^) 

The energy-dependent 1 Sq potential is given by 

v (p',p;E) = v LR (p',p) + X + 7 E. (2.6) 

Here vlr represents the long range part of the xPT potential, which is computed up to NLO or NNLO. However, 
the derivation presented below holds for any energy-independent long-range potential which is a function of p and p' 
and satisfies VL R (p',p) = VL R (p,p'). From now on we drop the partial-wave index, since this section only describes 
the singlet channel 1 So- 

To further simplify the presentation, we adopt the following operator notation of the LS equation 

t(E) = A + -fE + v LR + [A + jE + v LR ] g (E) t{E). (2.7) 
Setting E = in Eq. ([2~7]) leads to 

i(0) = A + v LR + [A + v LR ] a,(0) t(0). (2.8) 

Equation (|2.8|) contains only one unknown, A. Therefore, the matrix clement t(p',p;0) can be obtained from one 
experimental poi nt, i.e., the NN scattering length, a s , using the procedure explained in Subsection A, and given in 
detail in Ref. [331. Using properties of the LSE, Eq. (|2.8j) can be rewritten as 



t(0) = X + v LR + t{Q)g Q {0){X + v LR ) 

= [l + t(0). 9o (0)] [X + v LR ]. (2.9) 

This allows to reconstruct A if we wish to do so: 

A + VLB = [1 + i(0) 5 o(0)] _1 t(0). (2.10) 
Leaving E finite, Eq. (|2.7[) can be written as 

t{E) = [\ + 1 E + v LR ) [1 + g {E)t(E)} , (2.11) 
leading to a formal expression for the potential term 

X + 1 E + v LR = t{E) [1 + g (E)t(E)]- 1 . (2.12) 
Subtracting Eq. ((2~T0)) from Ea. ([2T2)) isolates the term -yE, 

1 E = t(E) [1 + goiE^E)}- 1 - [1 + <(0). 9o (0)]^ 1 t(0). (2.13) 

Proceeding in a similar fashion as in Ref. [Hj], we can obtain t(E*) at a fixed energy E* from the value of the NN 
phase shift at E* (for details see Appendix [A"|) . We find 3 

A + 7_E* + v LR = t(E*) [1 + goiE^tiE*)}- 1 = [1 + t(E*)g (E*)]-H(E*). (2.14) 



3 Note that some care is required if the standard subtraction method is used to deal with the singularity when the integration momentum 
becomes equal to po ■ ln that case it is useful to introduce an additional term in the integral equation that involves the potential evaluated 
at the point p' = po, p = po- However, if other techniques are used (e.g. contour rotation) this additional term is not necessary and so 
we do not include it in our derivation. 
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Taking the difference of Eqs. (f2~14|) and (j2~T2)) leads to 

1 (E* -E) = [l + tiE^goiE*)}- 1 t(E*) - t(E) [1 + g^E^E)]- 1 . (2.15) 
To eliminate 7, we insert Eq. (|2.13j) into Eq. (|2.15|) , 

(t{E)[l+g (E)t{E)T l [l + t(0)g (0)}-H(0)) - l) 

= [l + t^Xgo^r 1 ^^*)-^^) [l + 5 o(^W^)] _1 . (2.16) 

Rearranging Eq. (|2 . 16[) and multiplying both sides with 1 + gQ{E)t{E) from the right and 1 + i(0)<?o(0) from the left, 
we arrive at 

t(E) + t(0)[g o (0) - go(E)]t(E) + |^{t(0) - [1 + *(0)$o(0)] at(E*)} g (E)t(E) 

where a = [1 + t(E*)gQ(E*)]~ 1 , can be calculated numerically. Since we already obtained the fully off-shell i(0) and 
t(E*), Eq. (|2.17p is an integral equation for t(E), which we can solve by standard methods. 

To summarize, we perform two subtractions to the LS equation to eliminate the two unknown constants A and 
7. The resulting equation requires as input the scattering length a s and the phase shift at one specific energy E* . 
The only restriction on E* is that it must be within the domain of validity of our theory. Hence one can test the 
consistency of the theory by examining the extent to which results depend upon the choice of E*. 

C. Momentum-dependent contact interaction 

The coefficient of the momentum-dependent contact interaction, denoted here by C 2 , is not straightforwardly related 
to any S-matrix element. Thus we cannot apply our subtraction procedure in the case of a momentum-dependent 
contact interaction. Hence we instead adopt a "mixed" procedure, which involves a single subtraction plus fitting of 
C 2 . 4 

This "mixed" procedure is carried out as follows: First we guess a value for the constant C 2 , and then insert it into 
the half-off-shcll and on-shcll LSE: 

2 r A 

t(p,0;0) = v LR (p,0)+X + C 2 p 2 --M / dp' [v L r(p,p') + X + C 2 {p 2 +p' 2 )] t(p',0;0), 

77 Jo 

(2.18) 

t(0,0;0) = vlr(0, 0) + A — —M [ dp' (v LR (0,p') + A + C 2 p' 2 ) t(p', 0;0). (2.19) 

Taking the difference of the two equations cancels the constant A and leads to 
t(p, 0;0) - v LR (p,0)-v LR (0,0) + ^ + C 2P 2 

2 r A 

M / dp' (v LR (p,p')~v LR (0,p') + C 2 p 2 ) t(p',0;0). (2.20) 

77 Jo 

Using the already-chosen value of C 2 , together with the experimental value of a s , we can solve for i(p', 0;0) from 
Eq. (|2.20p . Then there is a consistency condition that determines the value of A. This equation can be easily derived 
from Eq. |(2T9|) : 

A = j - vlr(0, 0) + |M / A dp' (v LR (0,p') + C 2 p' 2 ) t(p', 0; 0) 
1-§M J A dp' t(p', 0;0) 



4 Note that this is different from the subtractive procedure of Ref. [45[ | for dealing with the same contact interaction. There the Born 
approximation is assumed to be valid for large negative energies. 
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The above equation gives the value of A that is consistent with the experimental value of the scattering length and 
our choice of C2. It thus defines a relationship A = X(C2',a s ). (Note that — in spite of the form of Eq. (|2.2ip — the 
relationship is not linear, since t(p', 0; 0) is also affected by the choice of C2.) Therefore, when trying to determine A 
and C2 we only need to guess C2 and can then obtain A from Eqs. (|2.20[) and (|2.2ip . These two constants are then 
entered into the on-shell LS equation which is solved for the phase shifts. Finally C2 is adjusted to fit the desired 
observable. In Sec. IIVI we will examine the results obtained when C2 is adjusted to reproduce the 1 So effective range, 
ro, and those found when we enforce the requirement that the theory correctly predict the phase shift at a particular 
energy E* . 

III. SUBTRACTIVE RENORMALIZATION IN THE J=l TRIPLET CHANNEL 

Since the NN interaction is non-central, the S = 1 waves constitute a coupled-channel problem. Here we consider 
the 3 Si- 3 Di coupled partial waves, and we again consider three different contact terms: 




Here we write the contact terms explicitly in their matrix form, where the diagonal represents the direct channels 
3 Si- 3 Si and 3 Di- 3 Di and the off-diagonal the channels 3 Si- 3 Di and 3 Di- 3 Si. The unknown constants arc A, 7 (or 
C2) and Aj, and their value is different in each case. Case A is the leading-order contact interaction discussed in 
Refs. @, H^. Cases B and C include the structures that appear at NLO [0(P 2 )] in the standard chiral counting for 
short-distance operators. 

At this order in the chiral expansion the J=l coupled-channel problem acquires two additional contact interactions. 
One is of the form <j\ q <T2 'Q [5(j an d hence gives a non-zero matrix element for the 3 Si- 3 Di transition. The other can 
be written as either an energy-dependent or momentum-dependent piece of the 3 Si- 3 Si potential, although until now 
most works on xET have considered only the latter [3, IE, 0, S| , but see Ref . [36[ . For the energy-dependent contact 
interaction we will show below that three subtractions can be performed to eliminate the unknown LECs. This leaves 
us with an integral equation for the i-matrix at an arbitrary energy that takes as input three experimental quantities. 
For case B, due to the momentum-dependence in the contact term, we use a double-subtraction-plus-one-fit scheme 
to solve the LSE. The procedures for solving case B and case C are quite similar, as they both employ subtractions 
to eliminate A and At. 

A. Constant contact interaction 

This case can be solved with a single subtraction, as described in Ref. [32|. Once again, the key fact is that once 
at, the triplet scattering length, is known, the value of i(0,0;0) is fixed, and that, together with knowledge of the 
long-range potential, is sufficient to obtain the two non-zero elements of t(p, 0; 0). The procedure by which we obtain 
these two elements, too and £207 can be reconstructed from Appendix [B] by setting C2 = in the first part of the 
derivation presented there. Hermiticity of v then implies that ioo(P)0;0) = too(0,p;0) and i2o(p, 0;0) = £02(0, p;0), 
and the same set of manipulations can then be used to obtain tmip' ,p; 0). The operator equation (|2.3p . with E* = 0, 
then allows us to reconstruct t(E), albeit this time as a two-by-two matrix of operators, rather than the single operator 
t that is needed to describe the 1 Sq channel. 

B. Momentum-dependent central part 



In this case we will eliminate two constants from the integral equation: A, and the coefficient of the tensor short- 
distance interaction At- We will then adjust the coefficient C2 to reproduce one piece of 'experimental' information. 
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As usual, we begin by writing the on-shell and half-off-shcll partial-wave LS equations: 

ho(po,Po;E) t 02 (p , Pa ;E)\ = f A + C 2 (Po + pl) + v° L R (p ,po) XtPo + v° L 2 R (p ,Po) 
ho (po , Po ; E) t 22 (po , Po ; E) J \ X t pl + v 2 j° R (p , p ) vf R (p , p ) 

+ 2 M f A M P ' 2 f X + C 2 (p 2 o+p' 2 ) + vI q r (po,p') \tP ,2 + vl 2 R (po,p') 
+ - Jo P Pl-P l2 + ie\ hP 2 + vl° R (p ,p>) v 22 R ( Po ,p') 

ho(p',Po;E) t 02 (p',p ;E) 
ho(p',Po;E) t 2 2(p' ', Po; E) 



(3.1) 



t 00 (p, Po ; E) t 02 (p,p ; E)\ _ ( \ + C 2 (p 2 +pI)+ vf R {p,pa) X t p 2 + v° L 2 R (p,p ) 



t 20 (p,po;E) t 22 (p,p ;E) 



Xtp 2 + vf R (p,p ) 



--M 

7T 



dp' 



P 



Po ~ P + ze 



A + C 2 (p 2 + p' 2 ) + vf R (p, p') X t p' 2 + vf R (p, p') 
hp 2 + vf R ( Pl p') vf R {p,p') 



*oo(p' ',Po;E) t 02 (p',p ;E) 
ho(p',Po;E) t 22 (p',p ;E) 



(3.2) 



where p\jM = E. The subscript (superscript) in tviiyjiR) indicates the angular-momentum quantum number of the 
channels. To simplify notation, we write the 2x2 matrix as t. Subtracting Eq. (|3.2p from Eq. (|3.ip cancels A, 



t(p,Po;E) - t(p ,p ;E) 



-M 



dp 



C 2 (p 2 - Pi) + vl° R (p,po) ~ v° L a R {p ,p ) v° L 2 R {p,p ) - v° L 2 R {po,p ) 
At (p 2 - Pi ) + vf R (p,p )- vf R (po , po ) vf R (p,p )- vf R Oq , po ) 
a 



V 



Po ~ P' 2 + * e 



C 2 (p 2 - pI) + vl° R (p,p') vl° R (p ,p') vf R (p,p') - vf R (p ,p') 
X t (p 2 p 2 ) + vf R (p,p>) v? R (p ,p>) v 2 L 2 R (p, P ') - vf R (p ,p') 



t(p',p ;E). 
(3.3) 



Letting E — -> 0, i.e. po — ► 0, in Eq. (|3.3p leads to 

( 



fo 0- 0) tfO 0-0) - ( C2p2 + 0) ~ U °°« (0 ' 0) ° 



-— M 

7T 



dv , ( C2P 2 + v™ R (p,p') - vl° R (0, P >) vf R {p,p') - vf R (0, P ') \ , 

X t p 2 + vf R {p,p') v 22 R (p,p') J*IP.°»°J- 



(3.4) 



Here we used the threshold behavior of the partial- wave potential, v l L l R (p, k) ~ p l k l , to infer that, e.g. vf^ R (0 : p) = 0. 

Equation (|3.4[) shows that this feature of vlr has as consequence that t 02 {p, 0; 0) = i 2 o(0,p; 0) = and t 22 (p, 0; 0) = 
^2 2 (0,p; 0) = 0. Using Eq. (|3 .4(1 we can thus obtain the first part of the t- matrix equation that needs to be solved, 
namely 



t 00 (p,0;0)-jj = C 2 p' + vl» R (p,0)-vl» R (0,0) 



--M 

7T 



dp' 



{C 2 p 2 + v? R (p,p') - v° L ° R (0,p')) tootf, 0; 0) + (v° L 2 R (p,p') - vf R {0,p')) t 20 (p' , 0; 0) , (3.5) 



where we have used that ioo(0,0;0) = j±. We observe that At has been eliminated from this equation, but that we 
do not have a closed set of equations, since t 2 o(p' , 0; 0) appears on the right-hand side. 
To eliminate X t we set p = in Eq. (|3.3p 

t(0,p ;E) - t(p ,p a ;E) = 

-C2PI + v° L ° R (0,po) - vf R {p ,p ) v° l 2 r (0,po) - vf R (p ,Po) 
- X t pl - vf R (p ,Po) - vf R (Po , Po ) 



-M 



dp' 



Po — P 12 + * e 



-CaP? + <k(0,p') - v° L R (po,p') vf R (0,p') - v 02 R {p ,p>) 



->>tP 2 o-v 20 R (po,p') 



-v 22 r(po,p') 



t(p',p ;E) 

(3.6) 
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In order to use the threshold behavior of £20 we need to divide Eq. (|3.3[) by p 2 — p$ and Eq. (|3.6[) by p 2 and then add 
both equations. A quick examination shows that 7 and A t both drop out of the final result — but the final result is 
quite lengthy. Hence, to simplify the notation arising from these manipulations we define the following map 



T[0(p,p')] = — L-5 
P-Po 



O(p,p')-O(p ,p')), 



(3.7) 



where the operator O can represent either the t-matrix or a potential function. With this notation we obtain 



F[t(p,p ;E)] - F[t(0, p ;E)] 



v°lr(p>Po) 
v 2 lr{p,Po) 

2 



T 
T 



7T 



^°fl(0,Po) 
vf R (0,po) 

,,'2 



M / dp' 



Pq — p' 2 +i£ \J~ 



T 
T 

T 



vf R (p,Po) 
vl° R (P,P') 



-T 

- T 

T 
T 



vf R (0,Po) 
v? R (0,Po) 

vTrM 

vIVo,p') 



T 
T 



v° l 2 r(p,p') 
v 2 lr(p,p') 



T 
T 



vf R (0,p') 
v 2 l 2 r(0,p') 



t(p',p ;E) 
(3.8) 



The second row of the kernel of this equation simplifies considerably since 



n<mpo)} = 



v 2 lr(po,Po) 



Pi 



due to the threshold behavior of vl R . A similar result applies to v 22 R . Thus we also have 

ho(Po,Po;E) 



F[t 20 {Q,po;E)] 



Pi 



(3.9) 



(3.10) 



Now, letting E — » (i.e. po — > 0) in Eq. (|3.8p . and extracting the S —> D element of the resulting equation, we can 
also employ (again because of threshold behavior) 



IukF[vI°r(p,p')] = 



(3.11) 



together with the analogous result for vjf R , in order to obtain 



W:P,0;0) 



P 

--M 

7T 



= lim 

po-»0 



ho(Po,Po;E) 



pI 



vl° R (p,0) 



— lim 

Po^O 



v lr(po,Po) 



Pi 



dp' 



V v 2 L ° R ( P0 ,p') 



lim 

Po^O 



Pi 



v 2 lr(p,p') 



lim 

Pa^O 



Pi 



too{p',0;0) 
t 2 o(p',0;0) 



(3.12) 



(3.13) 



This is the key result of our derivation, since it provides a second integral equation that, together with Eq. (|3.5[) . 
is a closed set for the half-off-shell t-matrix at zero energy. These equations require as input the long-range potential 
vlr, a choice of C2, and the experimental observables a t and the value of 



lim 

Po^O 



ho(po,Pa;E) 



Pi 



(3.14) 



This limit, and all the limits of Vlr in Eq. (|3 . 1 3[) . are well defined due to the threshold behavior of the potential and 
the t-matrix. The limit (|3.14p is a specific case of a "generalized scattering length" , which is defined as 



am 

—rr = hm 

M po^o 



ti'i(Po,Po]E) 



P l +l 



Since uy\ is related to an on-shcll clement of the t-matrix it can be obtained from an experimental phase-shift analysis 
(see, e.g. Ref. [52|V 
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Equations (|3.5|) and (|3 . 1 3|) . together with an initial guess for the value of C 2 , now yield too(p, 0;0) and t2o(p, 0; 0) 
by standard methods for the solution of integral equations. The other elements of the zero-energy t, t 02 (p, 0;0) and 
*22(p, 0; 0), are equal to zero, as mentioned above. 

It is formally possible to eliminate C 2 from Eq. (|3.8p by a similar set of manipulations to the ones that led to 
Eq. (|3.13p . This yields a pair of coupled integral equations for too(p, 0; 0) and t 2 o{p,0;0) that is NN LEC free. 
However, these equations contain the quantity lim Po ^o ^F[too{0,Po', E)\ as part of the driving term. This is a limit of 
a half-off-shcll matrix element, and so cannot be obtained in a model- independent way from experimental data. 

Therefore we now proceed similarly to the case of the momentum-dependent interaction in the singlet channel. 
We guess a value for C 2 , and solve Eqs. (|3.5[) and (|3.13|) for the half-off-shell t at zero energy. We then insert these 
quantities into the consistency condition for the LEC A t : 

A _ t 2 o(p, 0; 0) - v? R (p, 0) + §M / A dp' {vf R {p,p')t w (pi, 0; 0) + vf R {p,tf)t w tf , 0; 0)) 

P 2 (1-|A/ J o A dp% o (p',0;0)) 

Once we have At, the consistency condition for A reads 

A = at/M - vT R (0, 0) + |M J A dp' {{C 2 p>z + vT R (0,p'))t oo (p' , 0; 0) + (W 2 + vl\(0,p'))t 20 (p> , 0; 0))] 

1-fM f Q A dp't oa (p>, 0;0) 

The derivation of these two constraints for the LECs is given in Appendix [Bj Eqs. (|3 . 1 5[) and (|3 . 16[) define the full 
potential for the coupled-channels problem for a given trial value of C 2 . This potential is then inserted into the 
LSE and phase shifts computed. C 2 is adjusted until the desired experimental datum is reproduced with sufficient 
accuracy. 

This method has the advantage that Eqs. (|3.16[) and (|3.15[) define functional relationships: A = \(C 2 ;a t 7 a 2 o), and 
A t = Xt(C 2 ;a t ,a 20 ). (The generalized scattering length a 2 o enters implicitly through its effect on the half-off-shell 
t-matrices appearing in the formulae for A and At-) Thus, we need to perform only one-parameter fits in order to 
obtain all three of the LECs that are pertinent in the S-wave J=l channel. 



C. Energy-dependent central part 

Here we could follow the method of Section Hi Bl and eliminate all three constants from the integral equation: A, the 
coefficient of the tensor short-distance interaction, At, and the coefficient of the energy-dependent part, 7. However, 
it is simpler to mimic the steps of the previous Subsection, and reconstruct the underlying contact interaction from 
knowledge of the phase shifts. The manner in which this is done is described in detail in Appendix [B] but can be 
summarized in the following steps: 

1. Use knowledge of at, the triplet scattering length, to eliminate A from the integral equation. Hence obtain 
t oo (p,0;0). 

2. Substitute the result into Eq. (|3.4p . with 7 = (since the energy-dependent piece of the contact interaction 
has no impact on t at E = 0) and rearrange the equation so as to obtain At from the on-shcll quantity 
lim po ^o ho(po,Po;E)/pl = a 20 /M. 

3. The value of A can then be obtained from the requirement that ioo(0, 0; 0) = at/M. 

4. Use the 3 Si and 3 Di phase shifts and the mixing parameter e± at a specific value of the energy, say E*, to 
reconstruct t(pQ,pQ] E*). Then, from a suitably modified (|3. 3|) . solve for t(p,pQ-,E). 

5. With this half-off-shell f at a finite energy in hand, 7 can be obtained from the version of Eq. (|3 . 2[) that is 
suitable for the energy-dependent contact interaction. The final formula for 7 is given in Eq. (|B10[) . 

It may appear that in this procedure one needs five pieces of input data to fix the constants (at, the generalized 
scattering length a 2 o, and the three phase shifts at energy E*). However, only one of the four matrix elements of t is 
used to extract 7. Hence, in reality only three experimental inputs: £oo(0, 0;0), £oo(PoiPoi E*), and the derivative of 
t2o(po,Po', E) with respect to E at a specific energy, are needed to solve this problem. 

We then obtain phase shifts for arbitrary E by by substituting the obtained values of A, At and 7 into the potential 
that appears in the LSE and solving that equation for the on-shell t-matrix. 
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IV. RESULTS: THE SINGLET S-WAVE 

In this section we present the results obtained in the -"-So channel when we employ the three different types of contact 
terms introduced in Sect. II. First, we only adopt the leading-order contact term, i.e. a constant. Then, we consider 
the energy-dependent contact term, i.e., vsr,q = A + jE. Finally, we use the more standard momentum-dependent 
contact term vsr,o = A + Ci{j> 2 + p' 2 )- For each of these cases, we examine the results found with the following three 
different forms of the long-range potential v^r. 

1. TPE computed using dimensional rcgularization up to 0(P 2 ) (denoted as DR NLO). 

2. TPE computed using dimensional regularization up to 0(P 3 ) (denoted as DR NNLO). 

3. TPE computed using spectral-function regularization up to 0(P 3 ) (denoted as SFR NNLO). 

For the case of spectral-function regularization, we adopt A = 800 MeV as intrinsic cutoff. The effect that varying 
the intrinsic cutoff from 600 — 800 MeV has on this potential is discussed in reference [48| . When combined with the 
momentum-dependent contact interaction, the 0(P 2 ) TPE yields the usual NLO calculation in xET, while the two 
0(P 3 ) forms of vlr both give calculations of NNLO accuracy. 

A. The lowest-order contact term 

A proper renormalization of these TPE potentials entails the presence of contact terms up to second order in |q|. 
However, it is still of interest to investigate what happens if only the lowest-order contact term is included in the 
potential (see also Refs. [36|,[46|). This enables us to see how the phase shift converges as the LSE cutoff A is increased. 
Then, when higher-derivative contact interactions are added to the potential, we can see how the cutoff dependence 
changes. 

The renormalization for the leading-order contact term is performed with one subtraction, with the scattering length 
a s = —23.7 fm [49( used as renormalization condition. This is the same as procedure as described for the x So channel 
in Ref . [32| , except that here we are using a long-range potential that is higher-order in the chiral power counting. 

The phase shifts % calculated with the DR NLO, DR NNLO and SFR NNLO potential are shown for different 
cutoffs A in Fig. [U The results show that a single, constant, contact interaction is sufficient to stabilize the 1 Sq phase 
shift with respect to the cutoff once A > 800 MeV for the SFR NNLO TPE and A > 1000 MeV for the DR NLO 
TPE. The DR NNLO potential requires A > 1200 MeV to become stable. 

The DR NLO potential exhibits a resonance-like behavior when A < 700 MeV, but the phase shift quickly converges 
after A exceeds 800 MeV. This indicates that — at least for this potential — a cutoff smaller than 700 MeV removes 
essential physics, and brings the phase shift rather close to that obtained with vlb, = ^ope in Ref. [32j . On the other 
hand, it is amusing to note that after <5( 1 So) converges with respect to the cutoff, the DR NLO result is actually closer 
to the PWA93 phase shifts [13, |49| than is the result from the other two potentials. 

This indicates that the pure (i.e. un-compensated by contact interaction) NNLO TPE produces too large a correc- 
tion, although the general trend of providing more repulsion is correct. We note that the phase shift calculated with 
the SFR NNLO potential stabilizes as a lower value of A, and that this A — > oo result is closer to the PWA93 result, 
compared to that obtained with the DR NNLO potential. 

In the upper panel of Fig. [5] we plot the values of the effective range r obtained from the above potentials. Since 
a s is fixed to the same value in all calculations, this figure shows the influence of A on the low-energy phase shifts. 
Comparing our results for the three different potentials to the "experimental" value ro = 2.7 fm (49|, we see that the 
SFR NNLO potential gives a value of ro that is closest to the experimental value. The results with the DR NLO 
and SFR NNLO TPE are stable with respect to A once A > 700 MeV. The DR NNLO TPE value eventually also 
stabilizes within small variations once A > 1200 MeV (not shown in Fig. [2]). This cutoff is rather large, reflecting the 
too-strong energy dependence that was already discussed in the previous paragraph. 

B. Second order: energy-dependent contact term 

Next we use a second-order contact term with the simplest energy-dependent form, namely vsr,o = A + "/E. This 
form arises because the contact operator NHdtNN' N can absorb the divergences in the two-pion exchange potential, 
with the difference absorbed into terms that are proportional to the free-nucleon equation of motion. (See, e.g. Ref. 0], 
where the amplitude is calculated in Born approximation, and so the contact terms required are only a function of the 
nucleon on-shell momentum pq.) These additional terms can be transformed into higher-order multi-body operators. 
Thus, up to 0{P 3 ) in x?T, the energy-dependent contact term is just as valid a representation of the short-distance 
physics as the more commonly chosen momentum-dependent one. Of course, the use of energy-dependent potentials 
engenders additional complications — especially when one wishes to embed the two-nucleon Hilbert space in a larger 
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Fock space. Nevertheless, we can choose this form for the short-distance operators to describe NN scattering, and 
it has the advantage that the two unknown constants can be solved by our double-subtraction method discussed in 
Sect. II. 

As input this method requires the value of a s and the phase shift at a specific energy E* . In principle, this energy 
can be chosen arbitrarily. In our calculation we choose Ti a & = 2E* = 2.8 MeV, i.e. at the energy where the 1 So phase 
shift has its turning point. We tested other choices E* < 10 MeV and found that they have only a minor effect on the 
results. We also tested a relatively high energy, Ti a f, > 100 MeV. Such a choice generates results which fit the phase 
shift at low energies (due to a s as input) and at the chosen higher energy (due to the use of S(E*) as input). However, 
we did not adopt this procedure, since in it the ultimate dependence of results on the choice of renormalization point 
(i.e. E*) reflects the the cutoff-dependence in the phase shift in a manner that we shall discuss below. 

In Fig. [3] the % phase shift calculated with the TPE potentials DR NLO, DR NNLO and SFR NNLO with the 
energy-dependent contact term is displayed for cutoffs A = 600 — 2000 MeV. A general feature is that the phase shift 
oscillates as a function of A. There are also strongly diverging phase shifts at A = 800 and 2000 MeV for the DR NLO 
potential because a resonance state is created by the contact term. (A similar behavior for the DR NLO potential 
was observed in Ref. (3(| for a momentum-dependent contact term with A = 2000 MeV.) For the DR NNLO (SFR 
NNLO) TPE, at A = 1000(2000) MeV the phase shift diverges violently from the previous converged value. 

In order to further understand the effect of the energy-dependent contact term, wc calculate the DR NNLO TPE 
phase shift at Ti a b = 10 and 100 MeV for A up to 19 GeV and display the results in Fig. The phase shifts have 
the same oscillatory behavior, with the period becoming longer as A becomes larger. Comparison with the results 
obtained by adopting only a constant contact term (second panel of Fig. [TJ where this phase shift converges with 
respect to A for A > 1200 MeV), suggests that the quasi-periodic behavior of S(Ti a i, = 10 & 100 MeV) is caused by 
the energy-dependent contact term. 

The oscillatory behavior of the phase shift at this energy is associated with resonances at specific A. An argand 
plot confirms that there is a resonance for the cutoff at which this phase shift diverges. Tracking the energy of the 
resonance as a function of A shows that it moves in a regular way, with the real and imaginary parts of Er decreasing 
in size as A increases, until eventually the resonance becomes a bound state and the pole in t moves towards an energy 
of — oo. In order to understand this better, we perform a two-potential analysis of the t-matrix resulting from this 
particular vsr- The details are given in Appendix [Cl We find that the appearance of poles is a consequence of the 
specific form of the contact interaction. When these poles are in the domain of validity of the theory the resulting 
phase shifts are in poor agreement with the Nijmcgcn PSA. Only at those cutoffs where the artificial poles correspond 
to deeply bound states or high-energy resonances do the phase shifts match the Nijmcgen analysis well. 

The effective range ro obtained from all three different potentials is about tq — 2.68fm and varies by < 2% over the 
range of cutoffs considered, with the variation being only ~ 0.5% if the long-range part is computed to 0(P 3 ). This 
is a consequence of our choice of a low-energy phase shift as second input. Comparing with the upper panel of Fig. [2J 
one can clearly see that ro receives significant correction from the 0(P 2 ) contact term, especially for the DR NLO 
potential. 

Thus, the energy-dependent contact term improves the overall agreement of the ^o phase shift, and in particular 
the effective range ro, with respect to the average value ro = 2.7 fm extracted from the different Nijmcgcn potentials. 
Once one tries to increase the cutoff A beyond 1000 MeV, this form of short-distance physics leads to problems if the 
DR TPE is used as long-range potential. However, adopting TPE computed using spectral-function rcgularization 
means a wider range of cutoffs can be employed before the same issues occur as for the DR TPE. 

It is worthwhile mentioning that there exists a particular regularization of two-pion exchange which avoids resonant 
behavior in the phase shifts. This involves using dimensional regularization for the NLO part of TPE, and spectral- 
function rcgularization for the NNLO part 5 . Then, even if an energy-dependent contact term is chosen for the 
short-range potential, only minor oscillations of the phase shifts with respect to A occur for A > 1 GeV. However, 
strictly speaking, a potential which is obtained using two different rcgularization schemes is not self-consistent. Thus, 
in the later discussion we will only present the results obtained from the three TPE potentials mentioned at the 
beginning of this section. 

C. Momentum-dependent contact term 

A momentum-dependent contact term associated with the TPE has the form A + C2CP 2 + p' 2 ), and can be solved 
by one subtraction plus the use of a one-parameter fit, as explained in Sect. II. The input chosen for the subtraction 



5 see |35| for the detail of this potential and the analysis in p-waves. 
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is again a s . In what follows wc will try to perform two different fits to fix the value of Gi. One uses only low-energy 
information, namely the effective range, r = 2.7 fm. The other fit fixes C2 by attempting to reproduce the phase 
shift at T lab = 200 McV. 

The resulting 1 So phase shift is shown in Fig. [5] for the DR NLO TPE. In this case we cannot reproduce rp. or 
the phase shift at T lab = 200 MeV, once A > 600 MeV. This pathology of the fit at NLO was discussed in Ref. J4(|, 
and is related to the Wigner bound that limits the impact of short-distance physics on the effective range — or more 
generally on the energy dependence of phase shifts [42, [50] ■ Therefore, for the DR NLO TPE we can only perform 
an overall best fit to the phase shift. Fig. O shows that, within this limitation, the phase shift quickly becomes cutoff 
independent once A > 900 MeV. 

The results for the DR NNLO TPE are shown in Fig. [S[ where we see that the two different fit procedures generate 
different results for the same A. This is especially visible at A = 500 and 1000 MeV, where a resonance-like behavior 
is present in the latter case when C2 is fitted to r$. For values of A not close to these problematic cutoffs the phase 
shift is almost independent of the renormalization point. 

This independence of renormalization point is even more apparent when the SFR NNLO TPE is adopted. As can 
be seen in Fig. [71 for A between 700 — 1800 MeV, the two different fitting procedures lead to almost the same 1 So 
phase shift. The same figure also suggests that A < 600 MeV is too small and cuts off too much of the potential to 
allow a good fit to the Nijmegcn analysis. Finally, the agreement between the two fits breaks down at A = 2000 McV. 
Thus, by adopting the SFR NNLO TPE, we achieve renormalization-point independence for a wider range of cutoffs, 
but this renormalization-point independence breaks down for higher A. 

Finally, we show in the lower panel of Fig. [5] the extracted effective range ro as function of A for the 1 So phase shift 
fitted to the Nijmegen value at = 200 McV. This allows to assess the sensitivity of the results to renormalization 
point (E* f» or E* « 200 MeV lab. energy). The DR NLO TPE has the least variation of r with A over the range 
of A considered here, but it needs to be remembered that the experimental value cannot be reproduced with a real 
value of Ci 0, [5l|. At order NNLO, the SFR TPE gives values of ro which vary less with A than do the ones from 
the DR TPE once A > 800 MeV, and are stable until A = 2000 McV. This can already be inferred from Fig. [7J 

In summary, the results obtained from the NNLO DR TPE with a momentum-dependent contact term have a 
peculiar behavior at A = 1000 MeV. Adopting the SFR potential increases the validity range of A as in the energy- 
dependent case. 

V. RESULTS IN THE J=l TRIPLET CHANNEL 

In the J=l triplet channel we again adopt the three different potentials, i.e., the DR TPE up to NLO, the DR TPE 
up to NNLO and the SFR TPE up to NNLO. To renormalize the potential in this channel we employ three different 
contact terms (cases A-C in Sec. IIII[) . For case B (energy- dependent) and C (momentum-dependent), the value of 
the generalized scattering length 

lim ^2o(po,Po;-5') _ O20 
po^o pi ~ M 

is required. We first adopt the value given in [52| for 020- Then we adjust CX20 to find the best fit. Since 0120 is not 
given from experiment, we list the values extracted from several so called "high-precision" potentials in Table [U 

A. The lowest-order contact term 

We first perform the renormalization with the constant contact term (A) in the 3 Si— 3 Si channel. The results 
obtained with this short-distance potential, and long-distance potentials DR NLO, DR NNLO and SFR NNLO are 
shown in Figs.lSlfTUl respectively. In the first two cases the 3 Si phase shifts reach cutoff independence as A approaches 
1000 McV, and then diverge once more at higher A. The 3 Di phase shifts exhibit a similar feature, although the 
pattern of convergence is not as obvious as in 3 Si. For the mixing angle ei, there is no obvious convergence with 
respect to A for the DR NLO and DR NNLO TPE. This can be explained by the fact that we only adopt a contact 
term in the 3 Si— 3 Si channel and fit to the scattering length at = 5.428 fm. This was sufficient to stabilize the J=l 
triplet channel in the case of the one-pion-exchange potential[24|, [28|, [32||. However, the NNLO TPE is singular and 
attractive in both eigen-channels and thus needs two contact interactions (or, equivalently, two subtractions) in order 
to stabilize its predictions with respect to A 46]. In contrast, the NLO TPE has a repulsive behavior as r — > in 
both ,1=1 cigcn-channel. This also leads to predictions that are unstable after one subtraction, albeit for the opposite 
reason than at NNLO. 
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The SFR NNLO potentials do not have the same singular attractive behavior at short distances, and so it is no 
coincidence that its phase shifts show the best overall convergence in 3 Si, 3 Di, and e\. The SFR strongly reduces the 
short-range attraction in the NNLO pieces of the TPE, and a lower-order contact term is then able to remove the A 
dependence of the phase shifts. Fig. [TT] shows that the phase shifts at Tj j, = 10 and 50 MeV obtained with the SFR 
NNLO TPE and a constant contact term become stable for A > 2500 MeV. 



B. Energy-dependent contact term 

The results of using an energy-dependent contact term to renormalize the three different long-range potentials are 
shown in Figs. gMB In Fig. HI we adopt a t = 5.428 fm, a 20 = 2.28 x lO" 10 MeV" 3 [H2|, and the Nijmegen PWA 
values of S( 3 Si), S( 3 Di) and ei at Ti a b = 10 MeV as our input to renormalize the DR NNLO TPE. This choice docs 
not produce an optimal description of ei with respect to the Nijmegen analysis, although it is pretty close. 

Actually, ei is a notoriously delicate observable, and the choice made for a 2 o affects it appreciably. Therefore, in 
the case of DR NNLO we further adjust a 2a to 2.25 x 10~ 10 MeV -3 , which yields the best fit of ei with respect to 
the Nijmegen values for A = 700 — 1200 MeV. This result is shown in Fig. [T3J Comparison with Fig. fTJ] shows that 
this 1.5% shift in a 2 Q leads to noticeable improvement in t\ over this cutoff range, with few other significant changes. 
The same value, a 2 o = 2.25 x 10~ 10 MeV -3 , also gives the best fit for DR NLO, and this a 2 o produces results for 
SFR NNLO which coincide quite well with the Nijmegen data. In addition, u 2 § = 2.25 x 10 -10 MeV -3 is still within 
a few per cent of the value obtained for this quantity from the potentials constructed in Ref. [53| . Such a difference 
is certainly within the range of expected N 3 LO corrections to our results. 

The DR NLO results obtained with 020 = 2.25 x 10~ 10 MeV^ 3 are shown in Fig. [TU We see there that the phase 
shifts oscillate for A < 800 MeV, and then converge for A = 900 - 1800 MeV. For A = 2000 MeV, the phase shifts 
diverge badly again. 

The DR NNLO results present several interesting features. First, there exists a range of cutoffs A ~ 800 — 900 
MeV where all three phase shifts agree with the Nijmegen analysis remarkably well 6 . However, with higher cutoffs 
the phase shifts gradually move away from the Nijmegen analysis. Second, for A = 2000(1400) MeV, <5( 3 Si)(<5( 3 -Di)) 
diverge and show a resonance-like behavior. 

In order to analyze this feature, we calculate the J=l triplet phase shifts for a fixed energy of Tj a t, = 10 and 50 MeV 
and a wide range of cutoffs A = 500 — 5500 MeV, and plot them in Figs. [16] and |T7l Both figures show oscillatory 
behavior of the phase shifts with respect to A. The S and D-waves have different points of divergence in this cycle, 
and (unsurprisingly) affect one another when either gets large. The oscillation pattern of the mixing angle ei is a 
combination of that appearing in the two phase shifts S( 3 Si) and 5( 3 Di). 

In Fig. [15] we show the results for the SFR NNLO potential with a 20 = 2.25 x 10~ 10 McV~ 3 . The phase shifts only 
show minor cutoff dependence for A < 1000 MeV and converge nicely for higher A. As far as the fit to the Nijmegen 
PWA is concerned, a decent agreement is observed for all A for <5( 3 <Si) and 5( 3 Di). For ei, the converged curves do 
not fit as well as those for A = 600 — 800 MeV. This is not a surprise, since we adopted a value of a 2 o which fits t\ in 
this region of cutoff. We found that at higher cutoffs, i.e., A ~ 2300 MeV, the SFR NNLO TPE presents a divergent 
pattern in the phase shifts. Similar to the singlet channel, SFR NNLO can defer, but cannot avoid, the resonance 
caused by the energy-dependent contact term. 

The effective range r obtained at NLO is r = 1.75 ± 0.01 fm for 900 MeV < A < 1800 MeV, with the error 
representing the extent of the cutoff dependence in that range. For A = 2000 MeV, r$ = 1.70fm, which is due 
to the resonance-like behavior in the phase shifts created by the energy-dependent contact term. For DR NNLO, 
ro = 1.75 ± 0.01 fm except for A = 1400 and 2000 MeV, which are just the cutoffs where the phase shifts also diverge, 
as shown in Figs. [TBI and[T71 The same behavior in r$ is present in the SFR NNLO, but for higher cutoffs. 

In summary, in the triplet channel the energy-dependent contact term creates a similar pattern as in the singlet 
channel. For the DR TPE (SFR NNLO), there is a highest cutoff A = 1200(2300) MeV with which the potential can 
be iterated in the LSE before resonances are generated at energies within the domain of validity of xET. 



It is interesting to point out that for A = 800 — 900 MeV, the energy-dependent contact term actually produces phase shifts which fit 
the Nijmegen analysis better than the momentum-dependent contact term, in both 1 Sq and 3 Si — 3 D\ channels. This may indicate 
that this on-shell-energy-dependent contact term represents the short-distance physics quite well in this narrow range of cutoffs. 
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C. Momentum-dependent contact term 

The results of using a momentum-dependent contact term to renormalize the DR NNLO TPE are shown in Fig. [TBI 
Here we adopt a t = 5.428 fm and o^o = 2.25 x 10 -10 MeV -3 as input and then perform a fit to the phase shift d( 3 Si) 
at Tiai, — 200 MeV. This pins down the three unknown LECs in the potential at this order. 

For phase shifts obtained by our two-subtractions-plus-onc-fitting method, once the trial value of C2 is imposed, 
the other two constants A and A f are fixed by the two pieces of experimental information that are our input: 0:20 
and at- Thus, the results shown in Fig. [18] satisfy the constraints at Ti a {, = associated with this information and 
reproduce the 3 Si phase shift at T; a b = 200 MeV. However, they do not provide the best possible fit to the phase 
shifts. In contradistinction to the previous section, the fit is not appreciably improved by variations of 0120 • Compared 
to the results obtained using the same contact terms in Refs. @, 0, @]j the overall fit to the Nijmegen phase shifts 
is noticeably worse. Furthermore, the results obtained by fixing 5( 3 Si) at Ti a b = 200 MeV do not show stability for 
5( 3 Di) and t\ with respect to A. This suggests that the TPE associated with a momentum-dependent contact term 
depends strongly on the choice of renormalization point once A > 1 GeV. 

We wish to see up to what value of A it is possible to ameliorate these difficulties by a different choice of input. 
Thus we remove the two constraints imposed by a t = 5.428 fm and «20 = 2.25 x 10~ 10 McV~ 3 and allow all three 
NN LECs to vary freely, so that we can obtain the best overall fit to the phase shifts. The results for DR NLO, DR 
NNLO and SFR NNLO are listed in Figs. [T9l and [20l respectively. The fits are a noticeable improvement over Fig. [TBI 
For all three potentials, starting from A = 700 MeV, the 3 Si and 3 Di phase shifts fit the Nijmegen analysis quite well 
and have only minor variation with respect to A up to 1000 MeV. A stronger cutoff dependence shows up only in the 
mixing parameter e\. There we observe a clear decrease in the A sensitivity of the result as the long-range potential 
is improved from DR NLO to DR NNLO and then again to SFR NNLO. We plot rp extracted from those best fits 
in Fig. [5T] Again, we see better agreement with the Nijmegen value tq = 1.833 fm [52| for DR NNLO than for DR 
NLO. For 700 < A < 1000 MeV, r varies by less than 3% for all three potentials. 

VI. CONCLUSION 

In this paper we developed subtractive renormalization schemes which can be applied to NN scattering in xET 
at NNLO and maximize the use of information from experiment. We showed how to manipulate LSEs containing 
the contact interactions of xET so as to eliminate constant, tensor, and energy-dependent contact operators from 
the equation in favor of pieces of experimental data. A by-product of our analysis is the result that the momentum- 
dependent contact operator cannot be eliminated in this way, and thus is not straightforwardly related to any NN 
S-matrix element. We can still, however, analyze the xET potential containing that operator, by using a mixture of 
subtractive renormalization and more standard fitting methods. 

The subtractive-renormalization method shows particular advantages in cases where the LECs undergo limit-cycle 
or limit-cycle-like behavior, and thus adjusting several of them to reproduce data represents a difficult fitting problem. 
It is also very useful when carrying out calculations at cutoffs A > 1 GeV, since fine-tuning of the unknown NN LECs 
occurs in that regime. Finally, the direct input of experimental information allows a straightforward determination of 
the impact of the renormalization point — i.e. the choice of which experimental information we input — on the resulting 
phase shifts. 

Looking first at the case of a constant contact interaction, in the 1 So channel our results confirm the finding of 
Ref. : when a constant contact term is adopted, the resulting 1 So phase shift converges and becomes independent 
of A in the limit A 00. This occurs once A > 1000(1200) for the DR NLO and SFR NNLO (DR NNLO) potentials. 
However, we found that in this calculation the contribution from the NNLO(Q 3 ) part of the TPE is larger than the 
NLO part. 

In contrast, in the coupled 3 Si — 3 Di channels the dimensionally regularized interactions cannot be renormalized 
with a single, constant, contact interaction. However, when the two-pion exchange is computed using spectral-function 
regularization A-independent results are again obtained at large A: presumably because the short-distance behavior 
of this potential is the same as that of one-pion exch ange, w hich is well-known to be renormalized by a single constant 
contact term in the 3 Si — 3 Si part of the potential [24 l28l. l32| . 

We also considered the contact terms that can occur up to 0(Q 2 ) in xET, examining both energy-dependent and 
momentum-dependent contact terms. For both the 1 So and 3 Si — 3 Di channels, the phase shifts obtained from 
dimensionally regularized two-pion exchange show oscillatory behavior as a function of A when energy-dependent 
contact terms are used. This is associated with the appearance, and subsequent movement with energy, of poles in 
the NN amplitude. The first cutoff at which such a pole appears in the domain of validity of xET is around A = 1000 
MeV. 
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We also found that cutoffs less than about 600 MeV sometimes produced resonant features in the phase shifts. We 
attribute this failure to the propensity of these cutoffs to eliminate too much of the central attraction in the TPE that 
is responsible for reproducing the energy dependence of the S-wave phase shifts extracted from experiment — especially 
in the 1 So- 

In between these two limits, i.e. for 600 MeV < A < 1 GeV, we found that the 1993 Nijmegen PWA phase shifts, 
in both the 1 So and 3 Si — 3 Di channels, could be quite well reproduced using a short-distance potential that includes 
an energy-dependent contact term, together with a long-distance potential calculated to 0(Q 3 ) in x^T. With an 
energy-dependent contact term good fits to the Nijmegen phase shifts can be obtained at cutoffs above 1 GeV, but 
only in fairly narrow windows of A, where the energy of the spurious bound state produced by the energy-dependent 
contact term is large compared to the energy scales of interest. 

Finally, if the spectral-function regularized TPE is adopted, then the range of the cutoff one can use before the 
resonance-like behavior occurs is greatly increased. 

Turning to the case of the momentum-dependent contact interaction: we found that in the case of this type of 
contact interaction there is a significant improvement with respect to the Nijmegen phase shift when one improves the 
long-range part of the potential from NLO to NNLO. Indeed, if the xET NLO potential is employed there are certain 
observables (e.g. the effective range) that cannot be reproduced if we restrict ourself to real values of the LECs in 
the NN potential. If we wish our NLO calculation to reproduce the effective range in the 1 So we are restricted to 
A < 600 MeV. This problem is absent at NNLO. 

We examined the renormalization-point dependence of our results, and found that in the 3 Si — 3 D! channel at 
A = 800 MeV the results from both the DR and SFR TPE up to NNLO are almost independent of the renormalization 
point. However, once A = 1000 MeV is reached, the results with the DR NNLO potential show appreciable dependence 
on the choice of experimental input. Once again, the SFR TPE is better behaved in this regard: there is only minor 
renormalization-point dependence for the SFR NNLO potential from A = 600 MeV until A = 1800 MeV, and the 
renormalization-point independence breaks down at A = 2000 MeV. 

The use of a momentum-dependent potential and the expressions for the NNLO TPE should allow us to confirm 
previous results for the triplet channel (e.g. those of Ref. [5J]). However, we found that if we impose strict relations 
between the three unknown LECs in order to reproduce particular pieces of experimental data (e.g. at and 0:20)1 the 
resulting fit to the Nijmegen phase shifts is in general worse than that obtained by direct fitting of data over a finite 
range of energies [E[. Only when all the three constants are allowed to vary independently do we get best fits, which are 
comparable in quality to the results in [El [46|. However, we want be emphasize that this fit can only be obtained for 
cutoffs A = 700-900 MeV for the DR NNLO TPE. (For the SFR NNLO TPE, it is possible to obtain a reasonable fit 
up to A = 2000 MeV.) There is thus an appreciable difference between using two subtractions that fix the zero-energy 
behavior of the amplitude and performing a best fit involving all three J = 1, 5=1 NN LECs. That difference is an 
indication that the triplet S-wave phase shifts obtained from chiral TPE depend on the renormalization point rather 
strongly, despite the fact that they are quite cutoff independent in this window of A. 

In summary, we have probed the renormalization-point and cutoff dependence of chiral potentials which include 
TPE when these potentials are iterated using the LSE. We find that energy-dependent contact terms — which might 
be considered useful as they allow one to evade certain theorems regarding short-distance potentials [U, [42j — tend 
to introduce unphysical resonances in the results once A > 1 GeV. These resonances are generically absent when a 
momentum-dependent contact term is employed. Since the xET power counting suggests that energy- and momentum 
dependent contact interactions are equivalent up to the order to which we work this implies that the power counting 
is only a useful guide for A < 1 GeV. 

Although the momentum-dependent contact interaction does not produce the dramatic failures to agree with 
experiment that occur in the energy-dependent case, it does have some difficulties. Chief among these is the failure 
of certain aspects of the NLO calculation. Serious doubt must exist regarding the convergence of an expansion where 
(at least for the channel) LO is in poor agreement with experiment, NLO cannot reproduce the effective range, 
and it is not until NNLO that anything like a reasonable picture of the phase shifts emerges. 

The less singular interactions obtained with spectral-function regularization avoid many of the difficulties inherent 
in the use of their dimensionally regularized counterparts. The cutoff dependence is almost always weaker, with the 
cutoff dependence of phase shifts found with the SFR NNLO potential and a constant contact term disappearing 
altogether for A > 1 GeV. However, t\ cannot be reproduced in this range of cutoffs, and the issue with the effective 
range at NLO persists. 

We conclude that it makes little sense to iterate the two-pion-exchange potentials obtained from chiral perturbation 
theory using the Lippmann-Schwingcr equation unless one restricts the momentum in the LSE integral to lie below 
- 1 GeV [H [H B3. In this domain, and for the case of a momentum-dependent contact term, our subtractive- 
renormalization method produces results which agree with previous analyses [E|, HH, HE US Efil ■ 
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APPENDIX A: GENERATING t{E*) FROM A GIVEN VALUE OF 8{E*) 
For any given energy E* the on-shell value of the t matrix is related to the phase shift S(E*) by 

. ( * * pM foo(E*) e^ E 'hmS(E*) 

t(p ,Po,E ) = — r^— r = — , (Al) 

-Mp* -Mp* 

where foo(E*) represents the scattering amplitude and p$ = \J ME* . Then we write the half-shell and on-shell LS 
equations as 

t(p*',p* ;E*) = X + 1 E*+v LR (p*', P * ) + lM / A dp* p* 2 [ X + ^f 2 * + IT^J^ ) *(P*.PS;f )» (A2) 



7T 



Po — P + te 



t(p* ,p* ;E*) = X + 1 E*+v LR (p*,p*) + -M t dp* p* 2 ( X + 1 f* + l L 2 R ^ P * ] ) 

71 Jo V Po ~ P + l£ 

Subtracting Eq. (|A3|) from Eq. (TAT)) leads to 

e iS(E') s i n S(E*) 



t(p*',p*;E*) = v LR (p*',p*)-v LR (p*,p*) + 



-Mp* 

2 nsr f A J„* „*2 ( VLR(p*',P*) ~ VLr(P q,P*) 

V P*o 2 ~ P* 2 + i£ 



M I dp*p* 2 ( ULHyiJ jr~ 2 ULH r ,FJ ) t{p\pl;E*). (A4) 



The above equation does not contain any unknown quantities and can solved with standard techniques. With 
t(p*' ,Pq\E*) at hand, we can carry out another subtraction to get t(p*',k*;E*), where k* represents any arbitrary 
momentum. 

t(p*', k*;E*) = A + 7 £* + v LR (p*', k*) + —M [ A dp* p* 2 f X + ^+v LR (p*',p<') \ 



t(p* Q ,k*;E*)=X + 1 E*+v LR (p*,k*) + -M I dp*p* 2 ( X + 1 T + l L 9 R{P * 0,P * ) ) t(p*,k*;E*). (A6) 



7T Jo \ Po ~P *' + ie 

The difference of Eq. (|A5[) and Eq. (|A6|) gives 



t(p*',k*;E*)-t(p* ,k*;E*) = v LR (p*',k*)-v LR (p* ,k*) 

—M f A dp* P * 2 ( VLR b"f)-Z™te>rt)t(jr ik -.,E*). 

7T J V p* --p* 2 +ie J 

(A7) 



Using t(p*„ k*; E*) = t(k* ,p^; E*) wc can calculate t(p*', k*;E*) by solving Eq. (|A7|) . 
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APPENDIX B: DETAILS OF THE SUBTRACTION METHOD IN THE J=l TRIPLET CHANNEL 

Let us start with the on-shell and half-off-shell LSE in this, two-channel, case. We will consider a momentum- 
dependent short-distance potential, so we have: 

a t /M 0\ / A + u°° R (0,0) 
) " I 



2 r A 

-M / dp' 
* Jo 



X + C 2 p' 2 + v° L ° R (0,p') X tP ' 2 + vf R (0,p>) \ ( t 00 (p>, 0; 0) 
) \t 20 {p', 0;0) 



(Bl) 



( tooip, 0; 0) \ ( \ + C 2P 2 + v° L ° R (p, 0) 
^t 20 (p,0;0) J { X t p 2 +v 20 R (p,0) 

- 2 af f A ^f x + c ^p 2 +p ,2 ) + v lr(p>p') W 2 + vl\(p, P >)\ (tuVAO) o 
" * Jo V W + vl° R (p,p>) vf R (p,p') ) \t 20 (p\0;0) 

Subtracting Eq. (|B1[) from Eq. (|B2[) gives 

f t oo (p,0;0)-a t /M 0\ ( C 2 p 2 + v™ R {p, 0) 0\ 
^ ho(p,0;Q) 0) { X t p 2 +v 2O R (p 7 0) 0) 

2 M C2p2 + V ^ p ^- V ^^ v °lr(p>p')- v lR^p') 

( MfA0;0) o 

X {t 20 (p', 0;0) 
The only unknown in Eq. (|B3[) is At, thus we can express 

fA 



(B2) 



(B3) 



Xt = *2o(p, 0; 0) - vl° R (p, 0) + |M £ (^(p.pQtoo^, 0; 0) + v 2 L 2 R (p,p')t 20 (p', 0; 0)) 

p 2 (l-|A/ / A ^oo(p',0;0)) 



Similarly, once At is calculated, from Eq. (|Blj) we can express A as 

A = at/M v Q L R (0, 0) + |M / A ^ [(CV 2 + ^(0,y))t oo (y, 0; 0) + (X tP ' 2 + ^(0,^))^^, 0; 0))] 

1-fM J o A ^'t O o(p',0;0) 

Thus, by using the two scattering length a t and a 2 o together with a trial value for C 2 , we can obtain the other two 
unknown constants A and A f . This determines the complete NN potential at this order, and it is then a straightforward 
matter to compute the phase shifts using the LSE. 

The LSE with the energy-dependent contact term (case B) satisfies exactly the same relation as Eq. (|3.8[) and 
Eq. ([3~T5|) . while Eq. dSHJ) is replaced by 

Mp,0;0) - jj = vl° R (p,0)-vl° R (0,0) 
2 " A 



7T 







M dp' (vT R (p,p')-v° L o R (0, P ')) t 00 (p',0;0) + (v 02 R (p,p')-v 02 R (0,p'))t 20 (p',0 ] 0) 



(B6) 

Equation (|B6j) . together with Eq. (|3.8p . can be used to solve for t(p, 0;0). Note that here we don't need to input 
a trial value for C 2 . Once we have foo(P:0;0) at hand, we can substitute it, together with the t 20 obtained from 
Eq. (|3.13p . back into Eq. (|3.4[) (which is unaffected by the issue of energy- vs. momentum-dependence and so is the 
same as Eq. (|B4[1 ) and obtain At. 

The value of A then can be obtained from the new version of Eq. (|B1[) , 



(B7) 



a t /M 0\_(X + v°U0,0) 0\2 f A ,( A + vf R {0,p>) X tP ' 2 + v° L 2 R (Q, P ') \f t 00 (p>, 0; 0) 
0j-{ Oj n J Q dP { J {t 20 (p',0;0) 
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which gives 

A = at/M v%> R (0, 0) + |M / A dp' [v$> R (0,p')tao(p', 0; 0) + (W 2 + vf R (Q,p'))t 20 (p' , 0; 0)] 

/ A dp'ioo(p',0;0) 

Finally, in order to solve for 7, one needs a third piece of information as input. A convenient choice is the phase 
shifts at an energy E. From S\, 3 -Di and t\ at a specific energy E we compute t(po,po; E). Then, from a modified 
Eq. (O, 

t{p,p ,E) t ^P°> E )-{\ t (pi-pt)+vl%(p,po)-v? R (po,p ) vf R (p, Po )-vf R (p ,po) 

+ ^ 7o ^ P§ - P' 2 + fe I W - Pg) + wf«(p,y) - vf R (p ,p') v? R (p,p>)-v? R (p ,p>)) 

we solve for t(p,po\ E). Once we have t(p,po; E) at hand, 7 can be obtained from a modified Eq. ()3.2j) : 

too(p,Po;E) t 02 (p,p ;E)\ ( X + jE + v L ° R (p,p Q ) \ t pl + vf R (p, p ) 
ho(p,Po;E) t 22 (p,p ]E) J \ X t p 2 +vf R (p,po) v 2 L 2 R (p,p ) 

ln,r f A ^ P' 2 { \ + yE + v° L R (p,p>) \ t p' 2 +vf R (p,p') 



i /"A 

-M / dp' 
' Jo 



p 2 - p' 2 + IS V W + p') «Sk(p. 
too(p',Pa;E) t 02 (p',p ;E) 

t2o(p',po;E) t 2 2(p',po;E), 



as: 



A - t 00 (p,p 0] E) + vl° R (p,p ) + & M J Q A dp'p' 2 )]too( P \ P o;g+^^"»^^)^°^- p ° ;g ) 

L 7T JO ^ ^ pfi—p' 2 +lS ' 



(B9) 



(B10) 



Note that here we have chosen to solve 7 by matching the I = Z' = part of the matrix element in both sides of 
Eq. (|B9j) . One also could have chosen to solve 7 by matching the IV = 02 or the IV = 22 part of the equation. This 
corresponds to the fact that we have three phase shifts in the coupled triplet, namely 3 Si, 3 Di and s\, however we 
only need one additional condition to determine the value of the third unknown 7. 



APPENDIX C: TWO POTENTIAL ANALYSIS OF THE TRANSITION MATRIX AND ITS RELATION 

TO BOUND STATES 

Here we want to consider the TPE potential with an energy dependent contact term, i.e., 

\ + v LR {p',p)+ 1 E=V 1 +V 2 . (CI) 

where V 2 should contain the energy dependence. The general derivation of solution of the LS equation based on the 
Gellman-Goldberger relation [|35[ leads to the following expression for the ^-matrix 

t = h + (1 + hg ) V 2 (1 + g t) . (C2) 

Eq. (|C2[) can be solved by iteration, and the first order 

t {1) =t 1 + (l+t 1 g )V 2 (l+g Q t 1 ) (C3) 

gives the well known DWBA expression. However, we do not want to follow along that line. Instead we will use 
Eq. (|C3[) to derive a series representation of t in terms of V 2 and easily calculable integrals. Explicitly the result up 
to first order in V 2 , which we denote by t^ 1 ', is given as 

t { '(P,P,E) = h(p,p ;E) + V 2 (E) + - / dp p — 2 + - / dp p — 2 

7T Jo E-^M+ie tt Jo E-^j+ie 

+ i 2 -? [ A d^ p* 2 d P " p »*<r>*mwi<py-,m. m 
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Since V 2 = jE does not depend on the momenta p and pi it can always be taken out of the integral. Therefore, we 
define 



T{p;E) = - / dp*p 



* 2 h(p,p*;E) 



(C5) 



and observe that we may always make the replacement t\Ga — > F and G$ti — > T if the operators in question appear 
next to a V 2 . Using this notation we can rewrite Eq. (|(J4|) as 

t (1) =h + V 2 + TV 2 + V 2 T + TV 2 T. (C6) 

In the same notation, the sum up to second order in V 2 , t^ 2 \ will be 

i (2) = ti + (l+t igo )V 2 + {l + t m )V 2 g t^ 
= h + V 2 + TV 2 + (1 + hg )V 2 g t^ 

= h + V 2 + TV 2 + V 2 T + V 2 g V 2 + V 2 g TV 2 + V 2 g Q V 2 T + V 2 g t igo V 2 r 
+hgo [V 2 T + V 2 g V 2 + V 2 g TV 2 + V 2 g V 2 T + V 2 g ti go V 2 r] . 



(C7) 



Let us further define 



f(E) = {-f / dp*dpp* 2 p 



h(p,P*;E) 



{E -£± + ie)(E-^ + ie) 



2 f A 
G{E) = - dp 

I" JO 



(C8) 



and note that, when evaluated between two V 2 s, gotigo may be replaced by f{E). Consequently v- 1 ' and can be 
simplified as 



- ti + (i + r)(T/ 2 )(i + r) 
t (2) = h + (i + r)(v 2 + v 2 [G + f}v 2 )(i + r), 

The next order, reads 

t {3) = h + (1 + r) [v 2 + v 2 [G + f]v 2 + v 2 [G + f] 2 v 2 2 ] (1 + r). 

Continuing, the N t/l order t^- N ' can be written as 



tW = t 1 + (i + r) 



AT 



(i + r). 



Letting TV — > 00, we obtain the formal sum 

i = tl + (l + r) 



(i + r). 



(C9) 
(CIO) 

(Cll) 
(C12) 



This tells us that t— matrix diverges when ^ = f(E) + G(E). We have performed such a calculation for the bound 
states discussed in Sec. II V Bl and found that Eq. (|C12|) indeed holds at the energies where they occur. 
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Q?20 


2.28 2.28 2.09 2.18 



TABLE I: Value of the generalized scattering length Q20 extracted from various potentials, in units of 10 10 MeV 
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FIG. 1: (Color online) The 1 S NN phase shift as a function of the laboratory kinetic energy for different cutoffs A ranging 
from 0.6 to 1.2 GeV. The long-range potentials employed are DR NLO, DR NNLO and SFR NNLO (as noted in the y-axis of 
each figure), together with a constant contact term. The results are obtained by one subtraction with ao — —23.7 fm as input. 
The values of the Nijmegen phase-shifts [49| are indicated by the open triangles. 
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FIG. 2: (Color online) The effective range ro [in fm] in the 1 So channel extracted from calculations with the DR NLO (black 
square), the DR NNLO (red circle), and the SFR NNLO (green diamond) TPE combined with: a constant contact term (upper 
panel) , and a constant plus a momentum-dependent contact term (lower panel) . In both cases ro is shown as a function of the 
cutoff A in the LSE. In the lower panel, the coefficient of the momentum-dependent contact term is adjusted to reproduce the 
Nijmegen value of the phase shift at Ti a b = 200 MeV. The thick solid band represents the range of ro obtained from Ref . [52| . 
Note that for the upper(lower) panel, the value of r at A = 700 MeV for the DR NLO (SFR NNLO) potential is -10.2(0.15) 
fm, which is not plotted in the figure. 



24 



5 ( ls o> [ de §] [deg] 




FIG. 3: (Color online) The 1 So NN phase shift as a function of the laboratory kinetic energy for different cutoffs A ranging 
from 0.6 to 2 GeV. The potentials employed (from top to bottom) are DR NLO, DR NNLO and SFR NNLO, together with 
an energy-dependent contact term. The results are obtained by two subtractions with ao — —23.7 fm and the phase shift at 
Tiab = 2.8 MeV as input. The values of the Nijmegen phase-shifts [4S| are indicated by the open triangles. 
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FIG. 4: The 1 So NN phase shift at Ti a t = 10 MeV (upper panel) and 100 MeV (lower panel) as a function of the cutoff ranging 
from 0.5-19 GeV. The results are obtained using the DR NNLO potential with an energy-dependent contact term via two 
subtractions. 




FIG. 5: (Color online) The 1 So NN phase shift as a function of the laboratory kinetic energy for different cutoffs A ranging from 
0.6 to 1 GeV. The DR NLO long-range potential is employed, together with a momentum-dependent contact term. The results 
are obtained by making one subtraction with ao = —23.7 fm as input, and then per forming a best fit to the overall phase shift 
as given by the Nijmegen analysis. The values of the Nijmcg en phase shifts [49J are indicated by the open triangles. 
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FIG. 6: (Color online) The NN phase shift as a function of the laboratory kinetic energy for different cutoffs A ranging 
from 0.5 to 2 GeV. The potential employed is the DR NNLO TPE with a momentum-dependent contact term. The results are 
obtained by making one subtraction with ao = —23.7 fm as input and then performing a fit to either the effective range ro = 2.7 
fm (solid line) or the Nijmegen 1 So phase shift at Ti a t = 200 MeV (dashed line). The values of the Nijmegen phase-shifts [itj 
are indicated by the open triangles. 
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FIG. 7: (Color online) The 1 So NN phase shift as a function of the laboratory kinetic energy for different cutoffs A ranging from 
0.6 to 2 GeV. The potential employed is SFR NNLO with the momentum-dependent contact term. The results are obtained 
by making one subtraction with ao — —23.7 fm as input and then performing either a fit to the effective range ro = 2.7 fm 
(solid line) or the Nijmegen 1 So phase shift at Ti a b = 200 MeV (dashed line). The values of the Nijmegen phase-shifts are 
indicated by the open triangles. 



28 



[deg] A=600-ioooMeV [deg] A= 1200-2000 Me V 




FIG. 8: (Color online) The J=l coupled NN phase shifts as a function of the laboratory kinetic energy for different cutoffs 
A ranging from 0.6 to 2 GeV. The potential employed is the DR NLO TPE with one constant contact term. The results are 
obtained via one subtraction with ao = 5.428 fm as input. The values of the Nijmegen phase-shifts [13] are indicated by the 
open triangles. 




FIG. 9: (Color online) The J=l coupled NN phase shifts as a function of the laboratory kinetic energy for different cutoffs A 
ranging from 0.6 to 2 GeV. The potential employed is DR NNLO with a constant contact term. The results are obtained via a 
single subtraction with ao = 5.428 fm as input. The values of the Nijmegen phase-shifts is indicated by the open triangles. 




FIG. 10: (Color online) The J=l coupled NN phase shifts as a function of the laboratory kinetic energy for different cutoffs A 
ranging from 0.6 to 2 GeV. Here the potential is SFR NNLO with a constant contact term. The results are obtained by one 
subtraction with ao = 5.428 fm as input. The values of the Nijmegen phase-shifts are indicated by the open triangles. 
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FIG. 11: (Color online) The J=l coupled NN phase shifts at Ti a b = 10 (left panel) and 50(right panel) MeV as a function of 
cutoff ranging from 1-8 GeV. The results are obtained with the SFR NNLO potential and a constant contact term. 
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FIG. 12: (Color online) The J=l coupled NN phase shifts as a function of the laboratory kinetic energy for different cutoffs A 
ranging from 0.6 to 2 GeV. The potential employed is DR NNLO with a linear energy dependence in the central part of the 
contact term. The results are obtained by three subtractions with ao = 5.428 fm, a.20 = 2.28 x 10 -10 MeV -3 and the phase 
shift at Tiab = 10 MeV as input. The values of the Nijmegen phase-shifts are indicated by the open triangles. 
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FIG. 13: (Color online) The J=l coupled NN phase shifts as a function of the laboratory kinetic energy for different cutoffs 
A ranging from 0.6 to 2 GeV. The potential employed is the DR NNLO with a linear energy dependence in the central part 
of the contact term. The results are obtained by three subtractions with ao = 5.428 fm, «20 = 2.25 x RT 10 MeV~ 3 and the 
phase shift at Ti a t = 10 MeV as input. The values of the Nijmegen phase-shifts are indicated by the open triangles. 
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FIG. 14: (Color online) The J=l coupled NN phase shifts as a function of the laboratory kinetic energy for different cutoffs 
A ranging from 0.6 to 2 GeV. The potential employed is DR NLO with a linear energy dependence in the central part of the 
contact term. The results are obtained by three subtractions with ao = 5.428 fm, ct2a = 2.25 x KT 10 MeV" 3 and the phase 
shift at Tiab = 10 MeV as input. The values of the Nijmegen phase-shifts are indicated by the open triangles. 
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FIG. 15: (Color online) The J=l coupled NN phase shifts as a function of the laboratory kinetic energy for different cutoffs A 
ranging from 0.6 to 2 GeV. The potential employed is SFR NNLO with a linear energy dependence in the central part of the 
contact term. The results are obtained by three subtractions with ao = 5.428 fm, a.20 = 2.25 x 10 -10 MeV -3 and the phase 
shift at Tiab = 10 MeV as input. The values of the Nijmegen phase-shifts are indicated by the open triangles. 




1000 2000 3000 4000 5000 



A [MeV] 



FIG. 16: (Color online) The J=l coupled NN phase shifts at Ti a i = 50 MeV as a function of cutoff ranging from 0.5-5.5 GeV. 
The results are obtained with the DR NNLO potential and a linear energy dependence in the central part of the contact term. 




FIG. 17: (Color online) The J=l coupled NN phase shifts at Ti a b = 10 MeV as a function of cutoff ranging from 0.5-5.5 GeV. 
The results are obtained with the DR NNLO potential and a linear energy dependence in the central part of the contact term. 
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FIG. 18: (Color online) The J=l coupled NN phase shifts as a function of the laboratory kinetic energy for different cutoffs A 
ranging from 0.6 to 2 GeV. The potential employed is DR NNLO with the momentum-dependent contact term. The results 
are obtained by two subtractions with ao = 5.428 fm and «20 = 2.25 x 10~ 10 MeV~ as input and then performing a fit to the 
3 Si Nijmegen phase shift at Ti a b = 200 MeV. The values of the Nijmegen phase-shifts [49l | are indicated by the open triangles. 




FIG. 19: (Color online) The best fit for the NN 3 Si- 3 Di phase shifts as a function of the laboratory kinetic energy for different 
cutoffs A ranging from 0.6 to 1 GeV. The potentials employed are the DR NLO (left panel) and the DR NNLO (right panel) 
with a momentum-dependent central part of the contact term. The values of the Nijmegen phase-shifts [4!| are indicated by 
the open triangles. 




FIG. 20: (Color online) The best fit for the NN 3 Si- 3 Di phase shifts as a function of the laboratory kinetic energy for different 
cutoffs A ranging from 0.6 to 1 GeV. The potentials employed are the SFR NNLO with a momentum-dependent central part 
of the contact term. The values of the Nijmegen phase-shifts [4^] are indicated by the open triangles. 
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FIG. 21: The effective range ro [in fm] in the 3 Si channel extracted from calculations with the DR NLO (black square), 
DR NNLO (red circle) and SFR NNLO (green diamond) TPE and a momentum-dependent central piece of the contact term. 
Here ro is shown as a function of the cutoff A in the LSE, and is extracted from a best fit of the phase shifts to the Nijmegen 
PWA93. The solid line represents r obtained from Ref.[li). Note that the value of r at A = 600 MeV for DR NNLO is -11.5 
fm, which is not plotted in the figure. 



